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We consider model of quasi-2D magnetoelectric material as XY model for spin system on 
a lattice with local multiferroic-like interaction of spin and electric polarization vectors. We 
calculate the contribution of magnetic (spin) vortex-antivortex pairs (which form electric 
dipoles) to the dielectric susceptibility of the system. We show that in approximation of 
non-interacting pairs at T —> Tbkt (Berezinskii-Kosterlitz-Thouless temperature) dielectric 
susceptibility diverges. 


I. INTRODUCTION 

Magnetic and electric phenomena are finding more and more technological applications now. 
Therefore, materials that combine properties of ferromagnets and ferroelectrics (called multifer- 
roics) and materials that have coupled magnetic and electric subsystems (called magnetoelectrics) 
are interesting both from theoretical point of view and from point of view of possible technolog¬ 
ical applications (for reviews se^®^. Particulary interesting is to consider magnetoelectric and 
multiferroic thin filmJ^EI that contain topological defect^^. 

In the present work we consider magnetoelectric material, for example type-II multiferroic (ac¬ 
cording to classification, given irP), where effective interaction of electric and magnetic subsystems 
leads to electric polarization (P) and magnetization (M) coupling. We consider an easy-plane 
thin film (film plane coincides with the easy plane) with quite strong interplane interaction, so 
we assume that P and M vectors lie in the easy plane and their distributions are identical in all 
layers. We describe the magnetic subsystem of multiferroic in terms of classical two dimensional 
XY model for spin system on the lattice. Below Tbkt existence of a single vortex is energetically 
unfavorable, so that all vortices are linked in vortex-antivortex pair^. Coupling between polar¬ 
ization and magnetization brings us to the fact that the vortices in the XY model possess electric 
charges such that an electric charge is proportional to the topological charge of the vortexP. The 
resulting dipole gas contributes to the dielectric susceptibility of the material. In this paper we 
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calculate this contribution of magnetic vortex-antivortex pairs to the dielectric susceptibility of the 
magnetoelectric material. 


II. TWO DIMENSIONAL XY MODEL 

For magnetic subsystem of multiferroic we use classical two-dimensional XY model on the square 
lattice. Let Mj = M{cos (pi, sinc/ij) be magnetic moment at i-th site. Then the Hamiltonian of the 
system can be written as 

H =-j'^MiMj =-JM^'^cos{p)i-(pj). (1) 

(*.i> (*J> 

Here {i,j) denotes summation over all nearest neighbors and we consider ferromagnetic case when 
exchange integral J > 0. Taking continuous limit <pi —>■ (p{x) we get 

H=^Psl d^x{V(P)\ ( 2 ) 

where ps = JM‘^ is spin-wave stiffness and constant term is neglected. 

Apart from the ground state solution (p = const, there exist metastable vortex solutions 
(p = n arctan + (po- It can be shown that energy of single vortex with winding number n 

logarithmically diverge^I^ = nn^ps In ^ {R is of order of the size of the system, a is the lattice 
spacing), therefore single vortices don’t appear in macroscopic systems. 

There also exist metastable configurations, which are superpositions of single-vortex solutions, 
for example, vortex-antivortex pair configuration cp = ra(arctan ^Z'x ~ arctan ) ■ Energy of 
such vortex-antivortex configuration is finite: 

= 27rn^/95ln—. (3) 

a 

Here Y = ^{X, - Y,)2 + {Y, - y,)2 is the distance between vortex and antivortex cores. 

We see that single vortices don’t appear in a macroscopic system because of logarithmic di¬ 
vergence of their energy, but system of linked vortices and antivortices has finite energy. At 
temperatures T < Tbkt when vortices and antivortices are linked into pairs, there exists some 
equilibrium concentration of these pairs. 

III. INCLUDING MAGNETOELECTRIC COUPLING IN XY MODEL 

Consider the difference between pure XY model and XY model with coupled electric and mag¬ 
netic subsystems. In case of cubic lattice symmetry, keeping only the lowest-order terms, the 
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energy density of magnetoelectric material in electric field E is given 

p2 

w = - - 7 P((MV)M - M(VM)) - EP + a ^ diMj ■ diMj, (4) 

i,j=x,y 

where Xe is the dielectric susceptibility in the absence of M, 7 is the coupling constant, and 
|M| = Mq = const is the saturation magnetization. Minimization of Q with respect to P gives 

P = 7 Xe((MV)M-M(VM)) + XeE. (5) 

Let M(x) = Mo(cos i?i(x), sin(/>(x)). Then polarization is given by 

P = XXeMi ( ) + XeE. ( 6 ) 

\ dx(t> ) 

Inserting ([^ to Q, we see that electric, magnetic, and magnetoelectric parts of energy combine to 

w = {aM^ - Ixel^M^) - XejMi {d^cl^Ey - dycf^Ex) - \xeE\ (7) 

Term ^Xe-E^ gives constant contribution to the total electric susceptibility and further won’t be 
considered. Energy of one layer is 

H = a {aM^ - ^XeT^^o) J (f‘r{V(t)f‘ - aXe7^o J cPr{dx(t>Ey - dy4>Ex). (8) 

Assume, for a moment, that E = 0. In this case expression for energy ^ is similar to energy 
of one layer in XY model Q with effective spin-wave stiffness 

Ps = {2aMl - Xel‘^M^)a. (9) 

Hence, if we include the magnetoelectric coupling in conventional XY model, then in the absence 
of an external electric field we obtain the same XY model with different interaction constant. For 
typical values of parameter^ (oMq ~ 10“^ erg/cm, 7 Mq ~ 10“® (erg/cm) Xe = 1 — 10): aMg 
is greater than y^XeAfg by several orders of magnitude; therefore, ps remains almost unchanged 
by magnetoelectric interaction ps ~ 2aMQa. 

Let us discuss how magnetoelectric coupling affects the magnetic vortices. From Q it follows 
that polarization of magnetic vortex is P = — nyXeAfgr/r^. As it was shown iiP this leads to 
appearance of electric charge in the vortex core: 

Qe = 27rn7XeMo, (10) 

where qe is a vortex charge per unit film thickness. 
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Further, consider what happens in non-zero electric field. The second term on the right-hand 
side of (j^ becomes non-zero and ground state is reached when Vcj) = ^^{Ey, —Ex). This config¬ 
uration is a spin wave with period and it is perpendicular to the field. We see that A —>• oo 
as E ^ 0, hence vortex pairs of finite size aren’t influenced by infinitesimal field. But even 
at T = 0 spin wave creates polarization P = —s--—^-E and gives contribution to susceptibility: 

Ps 

Xsmnwave = ■ However, we expect that vortex pairs give divergent at Tbkt contribution to 

^ ps 

susceptibility and overcome nearly constant contribution of spin waves. Further, we don’t consider 
spin-wave contribution to susceptibility. 

When we apply an external electric field, boundary effects become important. Consider single 


vortex, which has electric charge (10) due to magnetoelectric interaction. From the total elec¬ 
troneutrality it follows that the boundary of the sample also becomes charged; it acquires charge of 
the same magnitude but opposite sign compared to the vortex. In an external field, the boundary 
charge effectively shields the charge of the vortex, and the value of screening depends on the geom¬ 
etry of the sample. In this work we assume for simplicity that our sample is a disk with radius R. 
For such disk screening reduces effective charge of vortices exactly twice (see Appendix for further 
details). 

We know that at any temperature below Tbkt there exist a certain amount of thermally 
activated vortex-antivortex pairs. Since vortex carries a positive charge and antivortex carries a 
negative charge, the pair forms an electric dipole. In the next section we calculate the dielectric 
susceptibility of such dipole gas with a variable number of dipoles. 


IV. DIELECTRIC SUSCEPTIBILITY CALCULATION 

Consider a system of non-interacting vortex-antivortex quasi-2D dipole pairs of vortex lines 
(at temperature T < Tbkt), which exist in thin film with thickness h, in electric field. Here we 
assume that h R and, therefore, all lattice layers have almost identical distributions of M, 
P, and electric charge density. For simplicity we consider only lowest-energy topological defects 
with n = ±I topological charges. Electrostatic energy of such dipole in external electric field is 
— ^grE (where r is the distance between vortex and antivortex cores and q'g = Qe/"^, see Appendix 
for details). Distance r can vary and this gives the contribution of vortex-antivortex interaction 
energy i) with effective spin-wave stiffness ([^ to the total energy. Hence for the total energy per 
unit film thickness of such a dipole we hav^ 

Utotal = -^e^E + 2ql^ In - - ij,, 

a 


( 11 ) 
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where = -Kri^ps/a', p is minus dipole formation energy per unit film thickness (i.e., —ph is 
an energy of a dipole with charges on neighboring sites). We consider the case of low dipole 
concentration, therefore, \p\h shonld be sufficiently large compared to Tbkt, which holds for XY 
model with very good accuracj®. Let U (r) = —q'^r^ + 2q^ In Y Then the grand partition fnnction 
is (in units ks = 1', (3 = ^) 


Z(E,T) = ^exp 


nph 


^ exp ( -h 


di f fevent 
configurations 


C/(ri) + ... + C7(r„) 


I / ^exp(-/3/i[/(r)) 


( 12 ) 


Here we replaced summation over all vortex confignrations by integration: Qcm is the coordinate of 
center of mass of a dipole. Let X be a number of lattice sites in one layer, then ex.p{/3ph) f ^ = 
Nexp{(3ph) and 


Z(E,/3) = Io 2JNexpi^ph) 


(Pi 


exp{—/3hU{r)) = 


(13) 


= /o(2v'Xi(E,/3)), 

where Io{z) = Yl’^=o ^he modified Bessel function of the firstll^kind and Xi(E,/3) is 

a partition function of one dipole: 

Z,(E,/J) =Wexp(/3,./.) / '^.^p{-lihU(r)) = 

J a 


drfd^Q 


r\-2l3hq^ + l 


exp (^[3hqf,rE cos (p) 
(14) 


Formally, this integral diverges at large r. However, we should keep in mind that we are calculating 
linear response (i.e., we take the limit E —0), so that r-divergence is cut off at radius of the sample 


r = R. Therefore, if we first take derivative of (14) with respect to E and after that let E = 0, 
then we obtain the correct result. 


Now, using (13), (14), in approximation of weak field E ~ 0, we can calculate contribution of 


vortex-antivortex pairs to dielectric susceptibility of the system: 




' = 


Sh dE^ 


TTO^q'lh'^ Nexp{Pph)Ii{25) VvNa^q'l exp{f)ph/2)^fihq^ — 1 • h{25) 


E=o 


ASh {hql,/2 - T)6 • h{25) 


AS 


ihql/2-T)-Io{26) 

(15) 


where S is area of the system, <5 = Temperature dependence 


X^f\T) is sketched in Fig. 1. 


Let us analyze applicability of formula (15) and find it’s asymptotics. When T —)• Tbkt 


hq^/2 susceptibility behaves as x — \j'-Tbkt \ ' dependence (and formula (15)) works up to 
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FIG. 1: Temperature dependence of vortex contribution into dielectric susceptibility Xe’\T) at T < Tbkt- 


temperatures T, at which vortex-antivortex pairs gas can be considered as dilute, namely when the 
average distance between pairs d = a-\/fdhq^ — 1 exp(/3|/r|/i/2) is much greater than the average si 


size 


of a pail® iX) = a ^ , This condition is satisfied in temperature region, where \T — Tbkt\ 

Vd^9m-2 

Note thaP /1 ~-vr^p^/a, so from /i ~ —vrq^. Hence, ~ — 


2 a 1^1 FV Tbkt 

Therefore temperature region |T — Tbkt\ — \'^bkt 


Tbkt 


%m/2 


= 27r. 


) ~ 10 ^Tbkt is a narrow region 


near Tbkt- This means that formula (15) works for all T < TbkT: that are not very close to 
Tbkt- 

(v) 

Let us find an asymptotic behavior of Xe at low temperatures. In this approximation T <C 
Tbkt ~ (which is equivalent to 2) and also /3|p|/i !:§> 1 (hence, 5 <C 1). Using 

asymptotic form Iu{z) ~ (i)^ /T{v + 1) for the modified Bessel functioiP^l^ -we obtain in zeroth 
approximation on E, retaining the principal T-dependent term: 

_ _ 2„,2]\/r2 / ,,h\ „2^,2^,2T\/r2 


x<”> = 5>rexp(/3;,ft)^ = | = exp 


T ) Aa- 2x‘^XeM^ 


(16) 


V. CONCLUSIONS 


In conclusion, in this paper we investigated properties of magnetoelectric thin film with easy- 
plane type of symmetry and type-II multiferroic-like interaction between electric and magnetic sub¬ 
systems. Magnetic vortices in such magnetoelectric possess electric charges and vortex-antivortex 
pairs form electric dipoles. Such dipole pairs have finite energy, so at any temperature T < Tbkt 
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there exists a certain amount of thermally activated vortex-antivortex pairs. We calculated the 
vortex-antivortex pairs contribution in static dielectric susceptibility of the system at temperatures 
T < Tbkt in approximation of non-interacting dipoles. This approximation is valid at tempera¬ 
tures which are not extremely close to Tbkt (when \T — Tbkt\ 3> \Tbkt exp(—In the low 
temperature limit dielectric susceptibility (15) behaves as an activation exponential ( |16[ ), which is 
consistent with the fact that at T —?■ 0 number of vortex pairs is proportional to exp As 


Tbkt formula (15) gives diverging susceptibility. This reflects the process of vortex-antivortex 


pairs unbinding and the phase transition. However, we expect that close to Tbkt interaction of 
dipole pairs becomes important and, therefore, at Tbkt susceptibility stays finite. 
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Appendix 


In this Appendix we calculate energy of charged vortex in electric field. As it was explained in the 
main text, magnetic vortex core acquires electric charge (10) due to magnetoelectric coupling (Q. 
Since our sample is electrically neutral, it’s edge acquires negative charge of the same magnitude, 


which should be taken into account (see Fig. 2a). 


Let our sample be disk with radius R, which contains single vortex with vorticity n, placed at 


distance Xq from the center of the disk. Choose coordinate system as shown in Fig. 2b, electric 
^eldE = {E^,Ey). 

From (IM polarization of single vortex configuration is 


P = 7XeMo^ 


-dy4) \ ^ n-iXeM^ I cose 

dxfp I ^ I sin 6* 


(A.l) 


Since we calculate linear response, XeE term in ([^ was omitted. Then electric energy of vortex 
configuration per unit him thickness is 

/ p27T r27r 

PE = nqXe-LfQ / d9 dr (Excose+Eysine) = nxXeMQ / d9{Ex cos 9+EySm9)l{e). 
Jo Jo Jo 

(A.2) 

Since l{9) = — Xq sin^ 9 — Xq cos 9 is an even function of 9, then odd term Ey sin 9 vanishes. 

Thus, 










8 



f 

- ' ^ y y-y'yfy*// f 

y yyy^/r/>//f f 
yyyy-'-//ft f 
'''■'■'//fft 
''''■//ftf^\\\\\\ 

yy^^/’/f>ti\\\\\\ 
//////(ft 
/ / / / / ^ I t I 
/III 


j\\\VV- 

\ \NVV; 
\ \ \ V 



a) Distribution of polarization for single vortex. 



FIG. 2: (a) Distribution of polarization for single vortex. Vortex core is positively charged, edge of the 
sample is negatively charged, (b) Coordinate system. Vortex core is placed at (0,0); disk center is (—Vq, 0). 


Uei = 2n^XeMQEx J dO cos 9 B? — Xq sin^ 6 — Xq cos 9^ = —TmxXe^oExXQ. (A.3) 
Using (10) we see that U^i = —IqgExXQ or for general displacement of vortex core r = (XqjIo) 


Uei = -q'eEr, 


(A.4) 


where gg = qel2 is an effective electric charge of the vortex, which is useful for finding electric 
energy of the vortex in an external electric field. 


Since polarization (A.l) is linear in phase angle cj), the same result (A.4) holds for vortex- 
antivortex configuration, with r denoting a vector that connects cores of the antivortex and the 
vortex. 
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